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Abstract 

Considering recently developed Hard Thermal Loop perturbation theory that takes into account 
the effect of the variation of the external field through the fluctuations of a conserved quantity 
we calculate the temporal component of the Euclidian correlation function in the vector channel. 
The results are found to be in good agreement with the very recent results obtained within the 
quenched approximation of QCD and small values of the quark mass (~ 0.1T) on improved lattices 
of size f28 3 x iV r at (N T = 40, T = 1.22b), (N T = 48, T = 1.452c), and (N T = 16, T = 2.98T C ), 
where N T is the temporal extent of the lattice. This suggests that the results from lattice QCD 
and Hard Thermal Loop perturbation theory are in close proximity for a quantity associated with 
the conserved density fluctuation. 
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I. INTRODUCTION 



Dynamical properties of many particle system can generally be studied by employing an 
external perturbation, which disturbs the system only slightly from its equilibrium state, 
and thus measuring the spontaneous response /fluctuations of the system to this external 
perturbation. In general, the fluctuations are related to the correlation function through the 
symmetry of the system, which provides important inputs for quantitative calculations of 
complicated many-body system. Also, many of the properties of the deconfined strongly in- 
teracting matter are reflected in the structure of the correlation and the spectral functions [l| 
of the vector current. 

The static thermal dilepton rate describing the production of lepton pairs is related to 
the spectral function in the vector current [2l, y] . Within the Hard thermal loop perturbation 

n h n 

theory (HTLpt) the vector spectral function has been obtained |2|, |4], |5|, which is found to 
diverge due to its spatial part at the low energy regime. This is due to the fact that the 
HTL quark-photon vertex is inversely proportional to the photon energy and it sharply rises 
at zero photon energy. On the other hand, the fluctuations of conserved quantities, such 
as baryon number and electric charge, are considered to be a signal js, 7] for quark-gluon 
plasma (QGP) formation in heavy-ion experiments. These conserved density fluctuations are 
closely related to the temporal correlation function in the vector channel through derivatives 
of a thermodynamic quantity associated with the symmetry, known as the thermodynamic 
sum rule It is expected that the temporal part of the spectral function associated 
with the symmetry should be a finite quantity and would not encounter any such infrared 



9] 



in 



divergence unlike the spatial part at low energy. A very recent lattice calculation 
quenched approximation has obtained the temporal part of the Euclidian correlation function 
associated with the response of the conserved density fluctuations, which is found to be a 
finite quantity. In view of this we would like to compute the temporal correlation function in 
the vector channel from the quark number susceptibility associated with the quark number 
density fluctuations within the HTLpt 10] and to compare it with recent lattice data joj in 
quenched approximation. 

The temporal correlator follows from the QCD polarization diagram. To lowest order 
perturbation theory it is given by the one-loop diagram containing bare quark propagators. 
The HTL resummation technique provides a consistent perturbative expansion for gauge 
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theories at finite temperature by using HTL resummed propagators and vertices As in 
usual perturbation theory it is strictly applicable only in the weak coupling limit but takes 
into account all dynamical effects. Going beyond the lowest order perturbation theory in the 
case of the temporal correlator, we use HTL resummed quark propagators and quark-gluon 
vertices in the polarization diagram (see 



fO 



12|). The HTL resummed quark propagators 



correspond to static external quarks (valence quarks). Within this approximation no internal 
quark-loops appear. In this sense our approximation can be compared to results from 
quenched lattice QCD. The inclusion of dynamical quark- loops requires to consider higher- 
order diagrams within the HTL resummed perturbation theory in which HTL resummed 
gluon propagators will show up. However, we do not expect a significant change in the 
result as higher contributions give less than 5% corrections in the case of thermodynamic 
quantities such as pressure [13I . Of course, close to the transition temperature higher order 
effects may become important jl^ . 

The paper is organised as follows. In sec. II we briefly discuss some generalities on cor- 
relation functions, fluctuation and its response (susceptibility) associated with conserved 
charges. In sec. Ill we obtain the relation between the response of the density fluctuation 
of the conserved charge and the corresponding temporal part of the Euclidian correlation 



function in the vector current. Next we compute them in HTLpt 10( and compare with 
lattice data. Finally, we conclude in sec. VI. 

II. GENERALITIES 

In this section we summarize some of the basic relations and also describe in details their 
important features relevant as well as required for our purpose. 

A. Correlation Functions 

The two-point correlation function jl-3] of the vector current, J M = ?/>(r, x)r A1 ^(r, x) with 
three point function T M , is defined at fixed momentum p as 

G^(r,p) = J d 3 x (J M (r,x)jt(0,0)> ^ ^ (1) 
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where the Euclidian time r is restricted to the interval [0,(3 = 1/T]. The thermal two-point 
vector correlation function in coordinate space, G^t, x), can be written as 

GW(r, x) = < J„(t, x) Jt(0, 0)) = rf | J^L G^K, p) , (2) 

n=— oo ' 

where the Fourier transformed correlation function Gfj, u {w n , p) is given at the discrete Mat- 
subara modes, w n = 2irnT. The imaginary part of the momentum space correlator gives 
the spectral function a(u, p) as 

0"ir(w,p) 



G^K.p) = - / da; 
=> cth(w,p) = -Im G ff (w,p) , (3) 

7T 

where H = (00, m, V) denotes (temporal, spatial, vector). We have also introduced the 
vector spectral function as oy = cr o + ^a, where an is the sum over the three space-space 
components and o"oo is the time-time component of cr MJ/ . 

Using 02]) and ([3]) in (Op) the spectral representation of the thermal correlation functions 
at fixed momentum can be obtained (2)] as 

Gh ^ V) = L dUaH ^ V) sinh[^/2] • (4) 

We note that the Euclidian correlation function is usually restricted to vanishing three 
momentum, p = 0, in the analysis of lattice gauge theory and one can write GuijT) = 
G h (t,6). 

A finite temperature lattice gauge theory calculation is performed on lattices with fi- 
nite temporal extent iV r , which provides information on the Euclidian correlation function, 
Gh{tT), only for a discrete and finite set of Euclidian times r = k/(N T T), k — 1, ■ • • N T . 
The vector correlation function, GV(tT), had been computed [2I within the quenched ap- 
proximation of QCD 1 using non-perturbative improved clover fermions 
bilistic application based on the maximum entropy method (MEM) [l6] 
N T = 16 and spatial extent N a = 64. Then by inverting the integral in d3j), the spectral 



15| through a proba- 



1 In comparison to thermodynamic quantities some quantities like mesonic correlation and spectral functions 
due to their structures are still too exhaustive and expensive to calculate in full QCD with improved 
lattice action. Hence the quenched approximation is still very useful to understand the various features 
of correlation and spectral functions. 
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function 2 was reconstructed 



HQ 



in lattice QCD. The vector spectral functions above the 
deconfinement temperature (viz., T = 1.5T C and 3T C ) show an oscillatory behaviour com- 
pared to the free one. In the high energy regime, u/T > 4 the vector spectral function, 



ay(u,0) agreed with that of the HTLpt 



BE 



On the other hand, the lattice spectral 



functions and dilepton rates |3|] were found to fall off very fast and became vanishingly small 
for u/T < 4 due to the sharp cut-off used in the reconstruction. 

n 

In a very recent lattice analysis [9J the low energy behaviour (viz., downward slope) of 
the spatial and vector spectral functions has been improved substantially on lattices up to 
size 128 3 x 48 by changing the slope upward through a fit to the free plus Breit-Wigner 
(BW) spectral functions at lower energy limit. This upward slope in the vector spectral 
functions in lattice gauge theory 9] resembles up to some extent that of the HTLpt spectral 
function {2], 4] and dilepton rate a, 18] at low energy regime despite the infrared problem of 
HTLpt at vanishing energy. Nonetheless, the existence of the van Hove peaks in the vector 
spectral function in HTLpt 0, [l9-21] has not been realized yet in the improved lattice 
analysis 9|, probably due to the ansatz that the lattice spectral function is proportional 
to that of the free plus BW one in the low energy regime. The existence of van Hove 
peaks cannot yet be ruled out, which are general features of massless fermions m a 

relativistic plasma in the low energy regime. Also, the high energy behaviour of the spectral 
function agrees well with those of free and HTLpt results, respectively, and are in conformity 
with its earlier analysis on the lattice 3|. On the other hand, the temporal component of 
spectral and correlation functions associated with the symmetry of the system are finite and 
do not encounter any infrared problem unlike their spatial part. Now, it would be interesting 
to analyse the temporal correlation and spectral functions, associated with the symmetry, 
within HTLpt and compare them with those of lattice gauge theory j^] . 



B. Density Fluctuation and its Response 

Let O a be a Heisenberg operator where a may be associated with a degree of freedom 
in the system. In a static and uniform external field T a , the (induced) expectation value of 



2 We note that lattice technique cannot directly be used to obtain the spectral function due to its difficulty 
to perform an analytic continuation of (J3)) from imaginary time to real time. 
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the operator O a (0, x^) is written jgl, [22I] as 

<l> a = (O a (0,*))r= 1 Tr \^/ W] l - = v jd*x(O a (0,lt)) , (5) 

where translational invariance is assumed, V is the volume of the system and l-i ex is given 
by 

The (static) susceptibility x»o- is defined as the rate with which the expectation value 
changes in response to that external field, 



Xaa{T) 



t=o 



(3 j d 3 x ((9 Q (0,^)a(0,^)) , (7) 



where (O a (0, 5c)O a (0, 0)) is the two point correlation function with operators evaluated at 
equal times. There is no broken symmetry as 

(o a (o,it))\^ = (aM))|^ o = o. (8) 

III. QUARK NUMBER SUSCEPTIBILITY (QNS) AND TEMPORAL EUCLID- 
IAN CORRELATION FUNCTION: 

The QNS is a measure of the response of the quark number density with infinitesimal 
change in the quark chemical potential, \x + 5fi. Under such a situation the external field, 
J^a, in OH]) can be identified as the quark chemical potential and the operator O a as the 
temporal component (Jo) of the vector current, J CT (i, x) = ipT a ifj, where T a is in general a 
three point function. Then the QNS for a given quark flavour follows from ([7]) as 

d 4 x /j (0,x)J (0,0)\ = -limRe G&(0,p), (9) 

\ / p->-0 

where Gq is the retarded correlation function. To obtain in concise form, we have used 
the fluctuation-dissipation theorem given as 

2 

G 00 (uj, p) = ~ l _ /T ImGg,(o;, p), (10) 
and the Kramers-Kronig dispersion relation 

Reg«(.,p) = r ^&(".fl , 

J~oo to W -U) 



dp 




d 2 V 




dfx 


(1=0 


dp 2 


-J 



where lim^o ImG^aj, p) is proportional to 5(u) due to the quark number conservation [8|, 
221 ]. Also the number density for a given quark flavour can be written as 

1 Tr [Me-W-^] (AT) _ dV 



P 



(12) 



V Tr [e-m-^O] V ' dn ' 

with the quark number operator, Af = f Jo(t,x) d 3 x = f ip(x)r i/j(x)d 3 x, and V = 
is the pressure and Z is the partition function of a quark-antiquark gas. The quark number 
density vanishes if fi — > 0, i.e., there is no broken CP symmetry. Now, ([7]) or fl9]) indi- 
cates that the thermodynamic derivatives with respect to the external source are related 
to the temporal component of the static correlation function associated with the number 
conservation of the system. This relation in is known as the thermodynamic sum rule js]. 

Owing to the quark number conservation the temporal spectral function <r o( w )0) in 
becomes 



1 



a 00 {u;,0) = -ImG&faO) = -uS(co) Xg {T) . 
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(13) 



Using ([TBI in (j4j), it is straight forward to obtain the temporal correlation function as 



Goo(rT) = -T Xq (T), 



(14) 



which is proportional to the QNS Xg an d T, but independent of r. The QNS has been calcu 



lattice gauge theory |23|-|3Q|], perturbative QCD |3lj, Nambu- 



lated within the framework o 
Jona-Lasinio (NJL) Model [22 
Ads/CFT correspondence and Holographic QCD 34J, Renormalisation Group approach 



321 ]. Polyakov-Nambu- Jona-Lasinio (PNJL) Model 



33 
35 



two loop approximately self-consistent ^-derivable HTL resummation 
I2I . I37I ] . We note that in a resummed perturbation theory 



11 



and HTLpt [10 



381 ] . the higher order loops 



contribute to the lower order due to the fact that the loop expansion and the coupling 
expansion are not symmetric. So, unlike conventional perturbation theory 39[ one needs 
to take a proper measure in order to calculate a quantity in a given order of a s correctly 



using HTLpt. We further note that various HTL approaches [12|, |36|, |37j, |41|, |42j have been 
used for calculating LO thermodynamic quantities and QNS in the literatures, which led to 
different results within the same approximation. Recently, we have developed a resummed 
HTLpt [10| by employing a variation of an external probe that disturbs the system only 
slightly from its equilibrium positions. In this way the effect of higher order variations of 
the external source is taken into account, which are essential to get the LO quantities cor- 



rect. Once this is done the LO quantities in HTLpt 



lOj agree with the HTL approaches of 



Ref. 36|, |42[. Here we would like to exploit the LO QNS in HTLpt obtained in Ref. 10[ for 



our purpose. 



Because of the structure of the HTL propagator [H|, [38( the LO QNS in HTLpt contains 
quasiparticle (QP) contributions due to the poles of the HTL quark propagator and Landau 
damping (LD) contributions due to the space like part of the HTL quark propagator. The 
QNS can be decomposed as 

xf TL (T)=x? P + X 9 LD (15) 
The LO QNS in HTLpt due to QP is obtained as 

Xq P (T) = AN c Nff3 j -j^L [ n ( u+ ) (1 - n(u + )) + n(w_) (1 - n(w_)) 

- n(k) (1 - n(k))] , (16) 

and the LD part is obtained [u]] as 



X L q u {T) = 2N c N f (3 J—Jdu \-r^) M". k ) n ( u ) C 1 - > ( 17 ) 

-it 

where n(y) is the Fermi-Dirac distribution, u + corresponds to the energy of a quasiparticle 
having chirality to helicity ratio +1, u;_ is the energy of a mode called plasmino having 
chirality to helicity ratio —1, and /3± are the cut spectral functions of the HTL quark 
propagator. The QP part results in (1161) is identical to that of the 2- loop approximately 



self-consistent ^-derivable HTL resummation approach of Blaizot et al 36|. The LD part 



([17]) cannot be compared directly to the LD part of Ref. [36| as no closed expression is given 
there. However, numerical results of the both QNS agree very well. We also note that Jiang 
et al [33] used HTLpt but did not take into account properly the effect of the variation of the 
external field to the density fluctuation, which resulted in an overcounting in the LO QNS. 
Moreover, in their approach an ad hoc separation scale is required to distinguish between soft 
and hard momenta and the thermodynamic sum rule is violated. In the HTLpt approach in 
Ref. [12I the HTL N-point functions were used uniformly for all momenta scale, i.e., both 



36] 



soft and hard momenta, which resulted in an overcounting within the LO contribution 
The reason is that the HTL action is accurate only for soft momenta and for hard ones only 
in the vicinity of light cone. 
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FIG. 1: (Color online) The 2-flavour scaled QNS with that of free one as a function of 
T/Tq- The solid lines are for LO in HTLpt whereas the dashed lines are for LO (propor- 



tional to g 2 ) in pQCD 



39, 



40l | . The different choices of the renormalisation scale are 



Q = 2ttT (red), and 4irT (blue). The symbols represent the various lattice data 



24 



28(]. The 



violet triangles (with T c = 204 ± 2 MeV), brown crosses (with T c 196 ± 3 MeV), green triangles 
(with Tr = 191 ± 2 MeV) and purple squares (with T c = 185 ± 5 MeV) represent p4 lattice QCD 
26], The squares (cyan) and stars (saffron) are from asqtad lattice QCD data j^j. The 



data 



23, 



solid circles (purple) represent quenched QCD data 



28l | for T c = 0.49A-pi7. The quark mass ranges 



between (0.1 to 0.2)m s , where m s is the strange quark mass near its physical value. Note that 



further lowering the quark mass to its physical value seems to have a small effect 
MeV. The details of these lattice results are also summarised in Ref. [3C|. 



29(] for T > 200 



Now we display in Fig. [T] the 2-flavour 3 scaled QNS in LO with that of free g 



unction o 
36 



39 



temperature that shows significant improvement over pQCD results of order q 



401 ] . Moreover, it also shows the same trend as the available lattice results [24H30|. 



We note that the QNS has a very weak flavour dependence that enters through the temperature dependence 



of the strong coupling as a s (T) 



127T 



(33-2JV / )ln(Q3/Ag) 



2t where Q is the momentum scale and Tc = 0.49Aq. 
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though there is a large variation among the various lattice results within the improved lattice 



(asqtad and p4) actions 



25 



271 ] due to the higher order discretisation of the relevant operator 



associated with the thermodynamic derivatives. A detailed analysis on uncertainties of the 



ingredients in the lattice QCD calculations is presented in Refs. [26|, |29|. This calls for 
further investigation both on the analytic side by improving the HTL resummation schemes 
and on the lattice side by refining the various lattice ingredients. 
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FIG. 2: (Color online) The scaled temporal correlation function with T 3 for Nf = 1 (left panel) 
and Nf = 3 (right panel) at T = 1.45Tc for Q = 2ttT (red) and 4-7rT (blue) as a function of scaled 
Euclidian time, tT. The symbols represent the recent lattice data e| on lattices of size 128 3 x 48 
for quark mass 0.1T in quenched QCD. 

Recently, an improved lattice calculation (9] has been performed within the quenched 
approximation of QCD where the temporal correlation function is determined to better than 
1% accuracy. Using the LO HTLpt QNS in ( ITS"]) we now obtain the temporal correlation 
function in (114)) and compare with the recent lattice data [9|. In Fig. [2] the scaled temporal 
correlation function with T 3 is shown for Nf = 1 (left panel) and Nf = 3 (right panel) at 
T = 1.45T C . We first note that the correlation functions both in HTLpt and pQCD have 



weak flavour dependence due to the temperature dependent coupling, ot s as discussec 



The LO HTLpt result indicates an improvement over that of the pQCD one 



36 



Defore. 
I for 
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different choices of the renormalisation scale as shown in Fig. [2j Also, the HTLpt result 
shows a good agreement to that of recent lattice gauge theory calculation {9] performed on 
lattices up to size 128 3 x 48 in quenched approximation for a quark mass ~ 0.1T. We also 
note that unlike the dynamical spatial part of the correlation function in the vector channel 
the temporal part does not encounter any infrared problem in the low energy part as it 
is related to the static quantity through the thermodynamic sum rule associated with the 
corresponding symmetry, viz., the number conservation of the system. We also presented 
two extreme cases of HTLpt temporal correlation function at T = 1.2Tc in Fig. [3] and at 
T ~ 3Tc in Fig.Hl respectively, for two different flavours and compare with the corresponding 
preliminary lattice data 43[, which are also found to be in good agreement. Finally, we 
also note that even if one compares _d , att ice action (a Sq tad) data and recent 
quenched data 9| for QNS, the quantitative difference is within 5% in the temperature 
domain Tq <T < 3Tc- 
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or online) Same as Fig. [2] but at T = X.ITq and the corresponding lattice data are 
4J with lattice size 148 3 x 40. 
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FIG. 4: (Co 
preliminary 



or online) Same as Fig. [2] but at T = 2.98Tc and the corresponding lattice data are 
4J with lattice size 148 3 x 16. 



IV. CONCLUSION 



The LO QNS as a response of the conserved density fluctuation in HTLpt when com 



pared with the available lattice data with improved lattice actions 24 



-2811 in the literature 



within their wide variation shows the same trend but deviates from those in certain extent. 
The same HTL QNS is used to compute the temporal part of the Euclidian correlation in 
the vector current which agrees quite well with that of improved lattice gauge theory cal- 
culations jol, l^] recently performed within the quenched approximation on lattices up to 
size 128 3 x 48 for a quark mass ~ 0.1T. It is also interesting to note that the quantitative 

nn 

difference between the recent quenched approximation data [9|, |43[ and the full QCD data 
with improved (asqtad) lattice action 27| for QNS is within 5% in the temperature range 
T c < T < 3T C . Leaving aside the difference in ingredients in various lattice calculations, 
one can expect that the HTLpt and lattice calculations are in close proximity for quantities 
associated with the conserved density fluctuation. 
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